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Abstract 

Financial markets based on Levy processes are typically incomplete and option prices depend 
on risk attitudes of individual agents. In this context, the notion of utility indifference price 
has gained popularity in the academic circles. Although theoretically very appealing, this 
pricing method remains difficult to apply in practice, due to the high computational cost of 
solving the nonlinear partial integro-differential equation associated to the indifference price. In 
this work, we develop closed form approximations to exponential utility indifference prices in 
exponential Levy models. To this end, we first establish a new non-asymptotic approximation 
of the indifference price which extends earlier results on small risk aversion asymptotics of this 
quantity. Next, we use this formula to derive a closed-form approximation of the indifference 
price by treating the Levy model as a perturbation of the Black-Scholes model. This extends 
the methodology introduced in a recent paper for smooth linear functionals of Levy processes 
[7] to nonlinear and non-smooth functionals. Our closed formula represents the indifference 
price as the linear combination of the Black-Scholes price and correction terms which depend 
on the variance, skewness and kurtosis of the underlying Levy process, and the derivatives of 
the Black-Scholes price. As a by-product, we obtain a simple explicit formula for the spread 
between the buyer’s and the seller’s indifference price. This formula allows to quantify, in a 
model-independent fashion, how sensitive a given product is to jump risk in the limit of small 
jump size. 
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1 Introduction 


The celebrated Black-Scholes model, which uses the geometric Brownian motion to describe the 
dynamics of the assets, is a cornerstone of the modern mathematical finance. However, it fails 
to reproduce significant features of empirically observed stock returns and option prices, such 
as fat-tailed distribution and implied volatility smile. For this reason, various extensions of the 
Black-Scholes framework have been developed in the literature. One popular approach is to replace 
the geometric Brownian motion with the exponential of a Levy process. 

Levy processes allow to quantify market risk much more precisely, but the option pricing 
problem for such processes becomes more involved. Exponential Levy models typically correspond 
to incomplete financial markets, meaning that the agents will not necessarily agree on a unique price 
for a derivative product. Instead, the price at which a market agent will accept to buy or sell a 
given derivative will depend on his / her risk aversion and preferences. A commonly used pricing 
paradigm in this context is the indifference pricing approach im which states that a fair price p of 
a contingent claim H for a market agent with utility function U and initial wealth Vo is the one at 
which the agent is indifferent between entering and not entering the transaction: 
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where S denotes the stock price and the maximum is taken over a suitable set of admissible trading 
strategies d. 

In this paper, we focus more specifically on the exponential (constant absolute risk aversion) 
utility function U(x) = —e~ ax , where a > 0 is the risk aversion parameter. This leads to a more 
explicit form for the indifference price: 


P a = ~ log ■ 
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exp ^— ol / q T dtdSt+aH^j 
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Additionally, using the so-called minimal entropy martingale measure (MEMM) denoted by Q* 
(see equation ©). the exponential utility indifference price can be expressed through a single 
optimization problem: 
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a. 


- log minE* 
a 0 


exp 



d t dS t +aH | 


( 2 ) 


where E* stands for the expectation under the measure Q*. 

Nevertheless, computing the utility indifference price ([2]) of even a simple European option in 
an exponential Levy model boils down to solving a non-linear integro-differential equation (see 
e. g ., mm), which is a tough numerical problem. This makes this approach unsuitable in a 
production environment of a bank, where prices must usually be evaluated in real time. For this 
reason, asymptotic approximations for the indifference price in incomplete markets are of great 
importance. 

One approach is to study the asymptotics when the number of contingent claims or, equivalently, 
the risk aversion a, is small [22 EH [23] EH m E]- In particular, for the exponential utility function, 
it is known EHEEUi] that under mild assumptions, as a tends to 0, the indifference price p a 
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converges to E*[i?], the expectation of the option’s pay-off computed under the MEMM, and that 
the optimal strategy converges to the quadratic hedging strategy under the MEMM. 

However, approximating the indifference price by the expectation under the MEMM fails to take 
into account the nonlinear features of the price. For this reason, in [22] , the authors compute the 
first-order correction to the exponential utility indifference price, and show that it is proportional to 
the residual risk of the quadratic hedging strategy under the MEMM: 


p a = E* [H\ 


— minE* 
2 0 


0 t dS t -H + E*[H] 


-f- o(o), ol —y 0. 


(3) 


These results are obtained under assumptions which is not straightforward to check for stock price 
models with jumps (in particular, Assumption 2 in [22]). Similar results for general path-dependent 
claims on a Brownian filtration are obtained in [29] using Malliavin calculus techniques. Our aim 
in this paper is therefore to obtain precise approximations for the indifference prices of options in 
exponential Levy models. 

First, in Theorem [T] we establish an approximation for the indifference price of the following 
form: 


p = E*\H] + - minE* 

J 2 ■d 
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/ "DtdSt — H + E*[ff] 
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- Error (a). 


(4) 


Unlike previous studies, our formula is non-asymptotic, in the sense that we provide an explicit 
bound on the error which is valid for all values of a smaller than a certain positive constant rather 
than asymptotically as a —> 0. In addition, this formula is proven under assumptions which are 
relatively easy to check in exponential Levy models. The proof of © is based on an interplay 
between the primal and the dual formulation of the indifference pricing problem to obtain an upper 
and a lower bound on the price, and is inspired by similar approaches in m and US- 

Next, we use formula Q to develop a closed form approximation to the exponential utility 
indifference price in exponential Levy models by treating the Levy model as a perturbation of 
the Black-Scholes model (see Theorem [2j) . In view of our non-asymptotic representation for the 
indifference price, this boils down to approximating the expectation of the pay-off under the MEMM, 
approximating the residual risk of the quadratic hedging strategy, and controlling the error term in 
Q. To this end, we use the approach suggested in a recent paper [?], which consists in introducing a 
one-parameter family of Levy processes (AT A ) t >o, where A = 1 corresponds to the model of interest, 
and A = 0 corresponds to the Brownian motion. One can then expand the quantities of interest in A 
around A = 0. In [7] such expansions were developed for the expectation and the residual risk of 
the quadratic hedging strategy in the case of European options with smooth C°° pay-offs. In our 
paper we employ a different technique to prove these expansions for a class of non-smooth pay-offs 
including for example the European put option. It is important to note that in problems of this 
type, regularity of the pay-off is an essential property which can strongly influence the convergence 
rate (compare for example with ]T51). 

As a result, we obtain an approximate formula for the indifference price of a European option in 
an exponential Levy model as a linear combination of the Black-Scholes price and correction terms 
which depend on variance, skewness and kurtosis of the underlying Levy process as well as high-order 
derivatives and a simple integral functional of the Black-Scholes option price. Our method is based 
on an interpolation between a Levy process and a Brownian motion and works well when the Levy 
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process in question is “not too far” from the Brownian motion. In a numerical study we compare 
our approximate formula to the exact value obtained by solving the integro-differential equation in 
the Merton jump-diffusion model and show that the precision of the approximate formula is good 
for realistic parameter values. 

Our approximate formula can be seen as an extension to the nonlinear utility indifference price 
of the valuation methodology based on the expansion around a proxy model, which was developed 
by E. Gobet and collaborators in the diffusion setting (see e.g., jj) and in [7] in the setting of Levy 
processes. It is also related to the “expansion in the model space” technique for utility optimization 
recently discussed in [25] . 

An important by-product of our study is a simple explicit approximate formula for the spread 
between the buyer’s and the seller’s indifference price: 


Ps ~Pb 


a 


Risk aversion 




Levy model 



Jump risk sensitivity of the option 


where a 2 is the variance of the Levy process at time 1, m 3 and 7714 are the third and fourth moments 
of the Levy measure and Pbs and E BS denote the option price and the expectation computed in the 
Black-Scholes model with volatility a. In our asymptotic regime, therefore, the spread is decomposed 
into a product of three factors: the risk-aversion which characterizes the economic agent, a factor 
depending only on the properties of the Levy model and a factor depending only on the variance a 2 
of the price process and on the properties of the contingent claim whose price is being computed. 
The last factor therefore provides a model independent measure of the sensitivity of a given product 
to jump risk in the limit of small jumps occurring at high frequency. 

The remainder of the paper is structured as follows. In Section [2] we present the precise 
mathematical setup for exponential utility maximization and indifference pricing in exponential Levy 
models. In Section [3] we derive a non-asymptotic approximation for the utility indifference price, 
and in Section [4] this formula is used to develop an expansion for the price “in the neighborhood of 
the Black-Scholes model”. Section [5] presents a numerical study of the performance of our expansion, 
and Section [ 6 ] analyzes specifically the formula for the bid-ask spread. Finally, proofs of technical 
lemmas are relegated to the appendices. 

2 Mathematical framework 

Exponential Levy models 

Let A be a Levy process on a probability space (fi, T, P), and let (Xt)t >0 be the completed natural 
filtration of A'. We fix a time horizon T < 00 and consider a price process S defined for t € [0, T\ by 
St = So£(X)t where So > 0 is a constant and £ denotes the Doleans-Dade exponential defined by 

£(X ) t = e x ‘-3W? (1 + AA s )e" AX '. 

0<s<J:AX s ^O 

The interest rate is taken to be zero throughout the paper. We make the following standing 
assumption. 
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Assumption 1. The process X is not a.s. monotone and there exists S < 1 such that X satisfies 
|AA't| <6 a.s. for all t £ [0,T]. 

Remark 1. The lower bound on the jumps of X ensures in particular that St > 0 a.s. for all 
t £ [0, T]. The non-monotonicity ensures that there exists a probability measure Q equivalent to 
P under which S is a martingale, which guarantees absence of arbitrage in the model (see e.g., /3 
section 9.5j). The upper bound on the jumps is a technical assumption needed in the subsequent 
developments. It is possible to assume that AX t < K for K > 1, however for notational convenience 
we impose the same bound on the negative and the positive jumps. 

Denote by (cr 2 ,^, 7 ) the characteristic triplet of A' associated with the truncation function 
x 1 —^ 1 |®|< 1 - By the Levy-Khintchine formula, this means that the for all t £ [0, T], the characteristic 
function of X t is given by 

E[e luXt ] = eif(u) = iu 7 — —f f ( e lux — 1 — iuxl\ [X \ < i)n{dx). 

2 J R 

Utility indifference pricing 

Consider a bounded contingent claim H = h(S) £ Ft, where h : O([0, T]) —> K is a mapping defined 
on the space of cadlag trajectories. When the financial market is not complete, this claim cannot be 
perfectly replicated, and therefore the price at which an individual agent will accept to buy / sell 
the claim will depend on the agent’s attitude towards risk, which may be quantified by a utility 
function. 

In this paper, we focus on the exponential utility function defined by U(x) = — e~ ax where 
a £ (0, 00 ) is the risk aversion parameter. Define the set of admissible trading strategies 

0 = {■!? £ L(S) | 3L* with E[e _ai *] < 00 s.t. (tf • S) t >L*Wt£ [ 0 ,T] a.s.} 

where L(X) is the set of ^-predictable A-integrable R-valued processes, and ( 1 ? • X) t := J^D s dX s 
denotes the stochastic integral with respect to X. Other definitions of the set of admissible strategies 
have been suggested in the literature m, but the above one appears sufficient in the context of 
exponential Levy models and it is somewhat more elementary and easier to check than the ones in 

US¬ 
As mentioned in the introduction, the seller’s utility indifference price of the claim H is defined 
by the implicit relation |l]), which, in the case of the exponential utility function, yields the explicit 
formula 

H _ U murage E [exp (-a(i? • S) T +aH T )} 
a ^ minaee E [exp (— a(d ■ 5)t)] 

The buyer’s indifference price p(( is defined in a similar manner and satisfies p(f = —p~ H . In the 
sequel, we shall focus on the seller’s price and omit the indices H and s. 

In conclusion, to compute the utility indifference price, we a priori need to solve two optimization 
problems. As shown below, in the context of exponential Levy models, the denominator of ([5]) can 
be computed explicitly. 

Let 

£(u) = 71 / H-—-b / ( e ux — 1 — ux)u(dx). 

2 Jr 
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Under the assumption that X is not a.s. monotone and has bounded jumps, £(u) is well defined for 
all u £ R, is bounded from below and there exists ip* £ R such that £(—atp*) = inf U £(u) (see the 
proof of Theorem 1 in [33]). Let d € 0 and (pt = fit St-. By Ito formula it is easy to show that 

M t = g-“ So VsdX s -J* t(-a<ps)ds 

is a local martingale. In addition, it is positive and bounded from above by e~ aL - Tln U^(«) 5 h ence 
a true martingale. Therefore 

E[ e -“(’ ? 'S')T] _ E[ e _ “ Jo Vt dx t-in i(-a Vt )dt+fT £(-acp t )dtj 

> E[e _ “^» T VtdXt-in e(-aipt) dt j e Jg e(-atp*)dt _ e Jg i(-OLtp*)dt . 

On the other hand, the strategy •&$ = is admissible. Indeed, as X has bounded jumps, its 
exponential moments are finite, and by Theorem 25.18 from Sato [5T], Efe - ^ lnf o<t<-r-XV] < qq. 
Thus, taking the strategy i9* we get equality in the above inequality, which shows that this strategy 
is optimal. Then we define 

dQ* _ e ~ a ^ XT 

dP E^ } 

Note that by definition of ip *, the measure Q* does not depend on a. Using Theorem 33.1 from 
[51] , it can be shown that under Q*, X is a martingale Levy process with diffusion component 
volatility er and Levy measure u*, where v*(dx ) = e~ a(p X u(dx). In particular, it implies that S is 
a Q* martingale. The measure Q* is the MEMM for S (see M)- Using this measure, the utility 
indifference price writes: 

p = - log inf E*e-“ ((,? ' 5)T - ff) . (7) 

a » 


Quadratic hedging 

Quadratic (also called mean-variance) hedging consists in finding an initial capital and a hedging 
strategy which minimize the expected squared P&L (Profit and Loss), that is: 

min E[(c + (i? ■ S)t — H) 2 ], 
ceiMee' L J 

where 0' is a suitable class of admissible strategies. We refer to m , isa, 0 for more details on 
quadratic hedging and to [20] for the specific setting of exponential Levy models. 

We shall see that the exponential utility indifference price is closely related to quadratic hedging 
under the measure Q*. Since Q* is a martingale measure, we can define 

Q' = {z? £ L(S) : d ■ S is a square integrable martingale}, 

and the optimal strategy can be computed as 

dt = ^5 > with H t =E*[H\X t }. 
d{S , S) t 

However, it should be noted that the optimal initial capital c may not be interpreted as a price 
of the claim H since it is equal to the price of the hedging strategy only and does not take into 
account the unhedged residual risk. 
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3 An approximation for the indifference price 

The goal of this section is to obtain a non-asymptotic approximation for the exponential utility 
indifference price in terms of the quadratic residual risk (error) under Q* in the Levy model under 
consideration. The quadratic hedging strategy under Q* will be denoted by d. 

Since X is a Levy process, the Q*-martingale (H t )o<t<r Las the predictable representation 
property [T 8 ( paragraph III.4d] and can be written as 

a s dX c a + f f 7 s{z)J x {ds x dz), (8) 

Jo Jr 

where X c is the continuous martingale part of the process X under Q*, Jx is the compensated jump 
measure of the process X under Q*, at is a predictable process and 7 S (-) is a predictable random 
function. 


H t =E*[H}+ [ 
Jo 


Theorem 1. Assume that there exists a constant L with 25La < 1 such that 

\H-E*[H]\ < L a.s., (9) 

|<Tf | < L a.s. for all t £ [0, T], (10) 

| 7 t(~)| < L\z\ a.s. for all t £ [ 0 ,T] and all z £ supp;/. ( 11 ) 


Then there exists a constant C a sL < oo such that for every e £ (0,1] the seller’s indifference price 
of the claim H satisfies 



f r T \ 2 

p-E*[H] - |E* 

U d s dS s -(H-E*[H})\ 


< a 1+s C aSL E* 


sup 

0 <t<T 


V s dS s -{H t -E*[H}) 


The constant C a SL can be chosen as 

C aSL = Ce AaL V (1 - 2 a8L)~ 2 , 
where C < oo is a universal constant. 

Remark 2. The formula of the above theorem is a non-asymptotic approximation formula for 
the indifference price, which can be used to recover a variety of asymptotic results. For example, 
observing that C a sL is bounded as a 0, we recover the asymptotics for small risk aversion: 


O' 


p = E*[H] + -E* 


d s dS s -(H-E*[H ]) 


+ o(a), a —> 0. 


Proof. Under the assumptions of the Theorem, 

- = 1 vcTt +J R z'Yt(zMdz) 

* S t - cr 2 + f R z 2 v(dz) 

and therefore \S t -'d t \ < L a.s. for all t £ [0, T]. We assume without loss of generality that E*[iL] = 0. 


2+e 
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Applying first Lemma [l] below with pay-off H' = aH and bound L' = aL, taking the logarithm 
using the inequality log(l + x) < x and dividing by a , we get 


p <|e- 


( rT \ 2 



/■* 

2+e~ 

/ J s dS s - H 

+ Ca 1+E e 4aL E* 

sup 

/ dsdSs - H t 


IV° ) \ 


0 <t<T 

Jo 
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Similarly, applying Lemma [2] below yields 


•s. ^ TO* 


( r T \ 2 



f t 

2+£- 

/ d s dS s - H 

- Ca 1+e { 1 - 2aSL)~ 2 E* 

sup 

/ d s dS s - H t 


\Jo ) 


0 <t<T 

Jo 
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□ 


Lemma 1 (Upper Bound). Let L be a constant such that assumptions ([9]) —(111 are satisfied. Then, 
there exists a universal constant C such that for every e £ (0,1], 



< 1 + X E* 

/ r T \ 2 



pt 

2+£- 

e -f?0 s dS s +H 

/ d s dS s -H\ 

+ Ce 4 i E* 

sup 

/ d s dS s - H t 


L 

2 

|_V° ) \ 


0 <t<T 

Jo 
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inf E* 

i960 


Proof. Introduce the stopping time 


r = inf{f > 0 : 


•d s dS s - H t 


> 1} A T. 


Since by assumptions, taking into account that S < 1, 


'dedSe — H 


< 3L + 1, 


the strategy = $ t lt< T belongs to 0 and we get 


inf E* 
tfee 


-ff# s dS e +H 


< E* 


e -ff#sdS s +H 


We shall use a Taylor formula of the following form: for every e £ (0,1) and m < oo, 

J2 


e^^l+rr-l- — + C\x\ 2+e , \/x £ [— m, m] 


with 


C = 




Thus, for C Le = (3L+1)1 r^ +1 


inf E* 

i?G0 



< 1 + -E* 

( [ d s dS s - h) 

+ c Le w 

[ d s dS s - H 

2+e 

L 

2 

\J o / 


Jo 



< 1 + -E* 


’dedSf, — H 



r T - 

2 +e" 

+ C Le E* 

/ d s dS s - H 
Jo 



+ (3 L + I) 2 ( - + C Le (3L + lfi )Q *[t<T\. 















































Then, by Markov inequality 




rt 




p 

2+e" 

Q* [ T <T]= Q* 

sup 

0 <t<T 

/ tfsdSs - H t 
Jo 

> i 

< E* 

sup 

0 <t<T 

/ tisdSs - H t 
Jo 



so that 


inf E* 

#£0 


e -f^ s dS„+H} 


< 1 + 




2 ' 


+ (Cie + 


(3L + 1) 1 2 
2 


+ C Le (3L + l) 2+e ))E* 



ft 

2+e" 

sup 

/ d s dS s - H t 


0 <t<T 

Jo 



Now, it is clear that one can choose a universal constant C such that the statement of the Lemma 
holds true. □ 


Lemma 2 (Lower Bound). Let L be a constant with 2 L5 < 1 such that assumptions © — ( |11| ) are 
satisfied. Then, there exists a universal constant C such that for every e £ (0,1], 


' ( f T - Y 

C 


[ - 

2 +£- 

U d s dS s -H\ 

E* 

(1 - 2LS) 2 

sup 

0 <t<T 

/ flsdSs - H t 
Jo 



Proof. From the results of j3] , we havc0 

P= sup |E Q [if] - -ff(Q|Q*)} , 

QeEMM(Q*) 


where EMM (Q*) denotes the set of martingale measures, equivalent to Q* and H( Q|Q*) is defined 

by 


H{ Q|Q*) = E* 



whenever this quantity is finite and equals +oo otherwise. Therefore, for any random variable D > 0 
such that DQ* is a martingale measure, 


p>E* [DH t } - E *[D log D\. 


( 12 ) 


Let k = | — LS, introduce the stopping time 


t k = infjf > 0 : 



d s dS s - H t 


> k} A T 


and define 

D = l+ f * ti t dS t ~H TK . 

Jo 

1 This reference provides a duality result for the class of admissible strategies which are bounded from below, but 

it can easily be extended to our class 0 using the dominated convergence theorem and the local boundedness of S. 
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By construction, E[D] = 1 and 

\D - 1| < k + \V Tk -S Tk -AX Tk | + | A H Tk | < K + 2LS < ^ + LS < 1. 
Moreover, for a bounded strategy d, 


E* 

D [ $ t dS t 

= E* 

D f K d t dS t 

= E* 

f f didSt-Hr) 

[ ^ d t dS t 


Jo 


. Jo 


[\Jo ) 

Jo 


because d is the optimal quadratic hedging strategy. Therefore, DQ* is a martingale measure. It 
remains to compute the right-hand side of (12). For the first term, using the Cauchy-Schwarz 
inequality and an estimate for r K , we get 


E *[DH t ] = E*[DH Tk ] = E* 


- / d t dS f 



( [T \ 2 


/ [T _ \ 2 

> E* 

H t — / VtdSt 

-E* 

H t — / d t dS t lr K <r 


\ Jo J 


\ Jo J 



( f T V 


r T - 

2 +e" 

> E* 

( Ht — d t dS t J 

-E* 

H t - d t dS t 

Jo 








< T] - 


> E* 


( r T \ 2 

- — E* 


rt 

2+e- 

H t - dtdS t 

sup 

/ d s dS s - H t 


V J 

K e 

0 <t<T 

Jo 



The second term in (12) can be estimated using the following Taylor formula: for every e £ (0,1) 
and A £ (0,1), 


x log x + 1 — x < 


C x-1 ) 2 
2 


+ C\x — l| 2+e , 


Vx £ [1 — A, 1 -(- A] 


with 


A 1-e 

6(1 - A) 2 ' 


Then, for C LSe 


_ (1/2 +L8) 1 ~ e 
— 6(1/2—L5) 2 ’ 


E* [D log D] = E *[D log D + 1 — D] < ^E*[(D - l) 2 ] + C LSs E*[\D - l| 2+e ] 


( r * \ 2 " 


r T k _ 

2+e" 

( H Tr - / d t dS t 

+ C LSe . E* 

H Tk - / dt.dSt 


\ Jo J 


Jo 



( r T \ 2 

+ (C L 8s + 4) E * 


ft _ 

2+£- 

[Ht- 0 t dSt 

sup 

/ d s dS s - H t 


\ Jo J 

K e 

0 <t<T 

Jo 
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Adding up the estimates for the first and the second term of (121, and choosing the universal 


constant C appropriately, the proof of the Lemma is complete. 


□ 
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Example 1. Let us check the assumptions of Theorem^ for the European put option with pay-off 
H = (K — St) + ■ The process ( H t )o<t<T is given by: 

H t :=E *[H\P t ] = P(t,S t ) 

where P(t, S ) = E*[(A' — S£(X)j’~t) + ], o-nd under suitable regularity assumptions on the process X 
(see e.g., m Proposition 2]), we have the martingale representation 

<T t dX( + / 7 s (z)J(ds x dz) 

Jo J R 

with 

a t = dP g^ St ~ St and Tt(z) = P(t,S t -{l + z)) - P(t, S t -). 

By dominated convergence: 

a dP(t, S) 

dS 

On the other hand, for z G supp v, 

\P(t,S(l + z))-P(t,S)\ <E*[|^(X) r _ t |l Sf(x)T _ t(1+ , )A1 <^] < 

4 Indifference price asymptotics in the neighborhood of the 
Black-Scholes model 

Since, as we have seen, the computation of the indifference price can be carried out under the 
MEMM, in this section, to simplify notation we omit the star in E*. In other words, we simply 
assume that all the expectations are taken under the MEMM unless specified otherwise, and that X 
is a martingale Levy process with diffusion component volatility cr and Levy measure v. 

In liquid financial markets, jumps are typically small and in most cases the Black-Scholes model 
provides a correct “order of magnitude” approximation to option prices. Thus it seems reasonable, 
in these markets, to treat more complex stochastic models as perturbations of the Black-Scholes 
price and to compute correction terms to this reference value. Our goal in this section is therefore to 
find an explicit approximation to the indifference price ([T]) in the situation when the Levy process 
X is “close” to the Brownian motion. 

To quantify what it means to be close to the Brownian motion, and following a recent paper by 
Cerny, Denkl and Kallsen [7|, we artificially introduce a small parameter A € (0,1) into the model, 
by considering the family of stochastic processes 

X t x :=XX t/x3 , 0 <t<T. 

Note that our parameterization is slightly different from the one introduced in [7] because that 
paper considers Levy models built using ordinary exponential, whereas we use the Doleans-Dade 
exponential. As a result, our formulas are somewhat simpler than the ones of [?]. 

With this parameterization, X 1 = X and, as A j. 0, X x converges weakly in Skorokhod topology 
to the process 

(A' t ) t >o = (dIT t )t> 0 , 


= E*[S£(A‘)T-tls£(.Y) T ^ t <if] < K. 


H t =E*[H}+ [ 
Jo 
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where W is a standard Brownian motion and a 2 = a 2 + f R x 2 B(dx). We then define: 

S x = S 0 £{X x ). 

Similarly, it is easy to show that the process S x convergences to So£(dW) as A tends to 0. 

Let H x = h(S x ) and consider the corresponding indifference price 

Px = - log inf E ( 13 ) 

a t?eO 


The following theorem provides an approximation of p\ when A —> 0 for European pay-offs, that 
is, we assume that H = Ii(St) and H x = h(S^). In this theorem and below, we let PBs{t,S ) 
denote the Black-Scholes price of the corresponding option computed with volatility d defined by 
d 2 = a 2 + J x 2 v(dx), that is, 


P B s(t,S)= E 


h 



^-{T-t)+sw T - t y ; 


where W is a standard Brownian motion. When t = 0 we also write Pbs( 0, S) = Pbs(S) to shorten 
notation. 


Theorem 2. Assume that 


The pay-off function h is a bounded, almost everywhere differentiable, the derivative h! has 
finite variation on [0, oo) and there exists L < oo such that \xh'(x)\ < L almost everywhere. 


J* x"^ is (dx~) 

• Either a > 0 or there exists (f £ (0, 2) such that liminf - > 0. 

rj .0 r 

Then, as A —> 0, 

Px = 


Pbs(S 0 ) + + ^fs^(So) 

{e^P^o) + 18S*PW(S 0 ) + 9S 5 0 P^ s (So) + So 6 ^ 6 J(So)} 


A 2 m 2 T 2 


^ ( m 4 - ^ ) E 


m 


?BS 


o2 d 2 p BS {t,s t ) 

* ds 2 


dt 


■ o( A 2 


where m 3 = f R x 3 is(dx), to 4 = f R x 4 i/(dx) and E ss denotes the expectation computed in the Black- 
Scholes model with volatility d. 

Remark 3. It is easy to check that the pay-off function of the European put option h(x) = (K — x) + 
satisfies the first assumption of the Theorem. Moreover, this assumption implies that assumptions of 
Theorem^ 7] (by an argument similar to the one given in Example [I]A As for the second assumption, 
it is satisfied by most parametric Levy models used in practice, such as CGMY f 6 j (with Y > 0) and 
normal inverse Gaussian W- It is not satisfied by the variance gamma model ]26f . 

Proof. The proof is based on the non-asymptotic approximation formula of Theorem [T] applied to 
the process S x , which takes the form 


p x -E[H x ]-^E 

( [ T d x dS x ~(H x -E[H x })\ 

< CE 

sup 

/W (2 *-£[**]) 

z 

lv° J J 


0 <t<T 

Jo 
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where i? A = argmin^ E 


(/q T fitdSt — (H x — E[if A ])y 


. The following lemmas provide estimates of 


the linear part of the price E[I7 A ], the nonlinear part of the price E 


(/q T ^sdSg — (H x — E[iJ A ]))' 


and the residual term in the right-hand side. In these lemmas we suppose that the standing 
assumptions of the paper hold true. Note that the expansion of the linear part of the price does not 
require the pay-off function h to be regular, but in the other two lemmas, regularity is an essential 
assumption without which the convergence rates as A 4- 0 may be different. 

Lemma 3 (Estimation of the residual term). Let the assumptions of Theorem^ hold true, let 
M A = f* dsdS* — (if A — E[/i(S 1 ^)]) and define Mj S = sup 0<t<T |M A |. Then Mq > 2, as A —> 0 


E [(Mj) ? ] = O A 9 (log 


A 


Lemma 4 (Estimation of the nonlinear part of the price). Let the assumptions of Theorem^ hold 
true. Then, as A —> 0, 


E 


- HSt) 


A 2 


m% 


= — to 4 —it E 


?BS 


s , 


2 d 2 p BS (t,s t y 2 

dS 2 


dt 


,(A 2 ). 


In addition, for the European put option with pay-off function h^Sr) = ( K — St) + , 


E 


BS 


d 2 P BS {t,S t ) 

dS 2 


dt 


K 2 

2nd 2 


e 1 +" du 
s/1 — u 2 


l oe £& _ r P 

where d = s 5 2 —. 

o\J± 

Lemma 5 (Estimation of the linear part of the price). Assume that 
• The function h is measurable with polynomial growth. 


• Either a > 0 or there exists ft G (0, 2) such that liminf 

r4-0 


x 2 v t dx ) 

' r 2-P 


> 0. 


Then, as A —> 0, 


E[/i(S A )] = P BS (S 0 ) 


Ato 3 T c3 p(3) ( a \ 

a ^0 ^BS W0y 


A mfT 4 p(4), q . 

94 ^Q^BSy^o) 


+ 


A 2 m 2 3 T 2 
72 


[6S 3 0 P^ s (S 0 ) + 18S 4 0 P^ s (S 0 ) + 9S 3 P^(S 0 ) + S 6 0 P^ s (S 0 )} + o( A 2 ). 


□ 
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5 Numerical results 


In this section, we illustrate numerically the performance of the asymptotic formula of Theorem [2] 
assuming that the asset price is described by Merton’s jump-diffusion model under Q*. Strictly 
speaking, this model does not satisfy the standing assumptions of the paper because the (log-normal) 
jumps are not bounded from above. However, in the numerical implementation discussed below, 
the Levy measure is truncated to a bounded domain (which can be chosen sufficiently large so that 
further increase of the domain does not modify the price). 


Merton’s jump-diffusion model In this model the stock price is defined by St = So£(X) t 
where 

N t 

Xt = /it + <r\Vt + 'y — 1) 

*=l 

where W denotes standard Brownian motion, jump sizes (Yi) ~ A/”(7,<5 2 ) are i.i.d. random 
variables and (N t )t> o is an independent Poisson process with intensity X AI accounting for the number 
of jumps up to time t. The Levy measure of X therefore has a density given by 


v(x) 


A M 1 x >_1 (log(tt + l)-7) 2 

- ——e 2 s? 

S(x + l)y/2n 


Implementation of the asymptotic formula In the numerical examples, we let A = 1 and 

approximate the indifference price by 


P =Pbs(So ) 


m 3 T q3p(3) (q x 


m 4 T 

24 


SZP ( b 4) s (S 0 ) 


m\T 2 
72 
aA 2 


■5 p(5) 
O r BS 


{e S 3 0 P™(S 0 ) + 18SqP B s(S 0 ) + 9 Si 

d 2 P BS (t,S t ) 


m 4 - ^ ) E 


BS 


s ? 


dS 2 


(So) 


dt 


c6 p(6) 

°O r BS 


(5o)} 


(14) 


Using the formula which has been justified asymptotically as A —> 0 for a finite nonzero value of A 
amounts to use a second-order Taylor expansion of a function at zero to approximate the value of 
this function at a point x ^ 0. The quality of the approximation does not depend on the specific 
value of x, but rather ou the smoothness of the function between 0 and x. The numerical examples 
of this section show that the indifference price is indeed smooth as function of A and that using the 
formula of Theorem [2] with A = 1 leads to a very precise approximation. 

To evaluate the approximate indifference price, one needs to perform three computations. 

• Evaluate a 2 and the moments of the Levy measure m3 and m 4 . In Merton’s model these 
quantities are easily computed from the explicit form of the Levy measure and are given by 

d 2 = a 2 + A M {e 27+252 - 2e^ + 1} 
m 3 = A M {e 3 ^ 52 - 3e 2 ^ +252 + 3e^ - 1} 
m 4 = A M {e 47+8<52 - 4e 37+ s 52 + 6e 27+252 - 4e 7+ T + 1} 
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Remark that although the original model has four parameters (since p is fixed by the martingale 
condition), the asymptotic formula only depends on three ‘group’ parameters a 2 , m 3 and 771 . 4 . 


Evaluate the integral in the last line of (14 ). In our example we consider the put option and 
evaluate the more explicit form of the integral given in Lemma [4] using an numerical integration 
algorithm. 


• Evaluate the derivatives of the Black-Scholes option price with respect to the underlying up to 
order 6. The exact explicit formulas for these derivatives are given in Appendix |D| 


Partial integro-differential equation and the finite difference scheme In this paragraph 
we briefly describe the HJB equation for the indifference price (see e.g., HUGH]) as well as the 
numerical scheme used to solve it. This scheme is inspired by well-studied schemes for linear integro- 
differential equations HU and is provided here only for the purpose of illustrating the asymptotic 
method; its full derivation and the study of its accuracy is out of scope of the present paper. 

Let Ht = (K — St) + and assume that S has the dynamics 


dSt^ 

S t - 


= dX t , 


where A is a martingale Levy process with Levy measure v and diffusion coefficient a. Then, the 
indifference price p(t , S ) satisfies (omitting the arguments where possible to save space) 


0 = 


dp S 2 c r 2 d 2 p 

— ^ ^ITds 2 


dt 


mm ■ 


+ J (p(t,S{ 1 + z))-p- Sz ^^j v (. dz ) 

( 1 ? - ||) +^J ( e <*(p(t,s(i+z))-p-sz#) _ 1 _ a (p( tj 5 (i + ,)) _ p _ Sz #?) v{dz ) | 


with terminal condition p(T , S) = (K — S) + . In log-variable x = log S , introducing P(t, x) = p(t, S ), 


dP a 2 (d 2 P dP\ 


dt 


2 \ dx 2 dx J 


dP 


0 =-ar + -S- «l 2 --flr +/ P(t,x + z)~ P ~(e z -I)-*- P(dz) 


dx 


rinj^ 2 + IJ (e a(p{t ’ x+z) - p - {eZ - lw ~l-a(P{t : x + z)-P-(e z ^l)d)^D{dz)\j, 


where v is the logarithmic transformation of v. 

To discretize this equation we introduce a time grid t, = ih , i = 0,..., N with h = j^, a space 
grid Xj = xq + jd, j = 0,, 2 M, and represent the Levy measure v as 

K 

u(dx) = ^ v k 5kd(dx), 

k— — K 

where K is an integer and 6 is the Dirac delta function. Let P l :l denote the approximation of 
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P(ti,Xj). We use the following implicit-explicit scheme: 


0 = 


P{ +1,7 Pi. 'i , ® ( Pi,i—1 T Pi ,7 + 1 2 Pi j Pi 7-1-1 Pi, j— 


d 2 


j,j+1 r uj-l 
2d 


+ J2 (Pi+i, j+ k-Pi + ij-(e kd - 1 )^ 

k——K ' 


j' + l \ - 


2 d 




+ min ■ 


0(7“ 

2 


d - 


Pi 


t+l,j+l 1 

2 d 


+ “ H ^e“^ p(i+1 ’2+ fe ) _p ( i + 1 d)-( ei!<i - 1 )' ? ) _ i _ a (p(i + l, j + fc) - P(j + l, j) - (e fcd - l)d)) 




k——K 

In other words, introducing the notation 

K 


k— — K 


B j (P i+ i)= [Pi+ij+k-Pi+ u-(e“-l) 


Pj+lj+l Pi+l,j —1 \ - 


2 d 


Z'fc 


and 


^■(P i+1 ,d) = — d- 


P: 


i-f- l,j + l P+1 ,j — 1 

2 d 


+ i 5Z (e a(p t i+1J+fe P p ( i+1 ’ J 'P( e '“ i - 1 ^) - 1 - a(P(i + 1, j + k) - P(i + 1 ,j) - (e kd - l)d)) D k , 


k——K 

we have for j = 1 ,..., 2M — 1 


p,t (i 


d 2 




oh o‘ 

2 d 2 + Id 


- P 


*,j + l 


(7/1 (J 

~2d? 4d~ 

= P+i j + hBj(P i+ 1 ) + /iminPj(P+i,d) 

$ 


with boundary conditions 

Pi, 2 M = (K — e X2M ) + and P, 0 = (tf - e x °)+. 

Numerical comparison 

In this paragraph we compare numerically the asymptotic formula for the indifference price of a 
European put option with the value obtained by solving the PIDE using a finite difference scheme. 
The computational time required to solve the PIDE with adequate precision is 107 seconds on an 
iMac with 2.8 GHz Intel Core i7 processor (the implementation was done with Python programming 
languate, using a single processor core). The parameters of the scheme were N = 40 (number of 
time steps), 2 K = 100 (number of points to discretize the Levy measure) and 2 M = 200 (number of 
space steps). 

Figure [T] left graph, plots the price computed with the two methods as function of the initial 
price of the underlying with risk aversion parameter value a = 10. For comparison, we also plot 
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Figure 1: Left: Indifference price computed with PIDE and with the asymptotic formula, in Merton 
model (log-normal jumps), as function of the initial price of the underlying So for the risk aversion 
parameter a = 10. Right: Indifference price computed with PIDE and with the asymptotic formula 
as function of the risk aversion parameter a for So = 1. Other model parameters: strike K = 1, 
maturity T = 1 year, diffusion volatility a = 0.2, jump intensity X M = 5, average log jump size 
—5%, log jump size standard deviation 10%. 


the linear part of the price (E*[iJ]), computed using the explicit formula available in the Merton 
model. As can be seen from the graph, the bid-ask spread (that is, twice the difference between 
the indifference price and the linear part of the price) for an at the money option corresponds to 
about 6% of the option price. This is a rather high value for the spread, which means that the risk 
aversion parameter value which we use is also rather high. 

The graph clearly shows that for the chosen parameter values, which correspond to a realistic 
market scenario, the asymptotic formula is quite precise. To further explore the domain of validity 
of the approximation, in the right graph of Figure |T] we plot the price of an at the money put option 
(that is, we take Sq = 1) as function of the risk aversion parameter a with a higher resolution. 
Remark that the asymptotic formula for the indifference price is linear in a. We see that in this 
example the asymptotic formula reproduces the linear component of the price with almost no error 
(for a = 0), and the nonlinear component of the price with high precision, even for relatively large 
risk aversion values. 


6 Bid-ask spread and sensitivity of options to jump risk 

As a by-product of the asymptotic formula of Theorem [2] we obtain a simple explicit approximation 
for the difference between the seller’s and the buyer’s indifference price of a European option, that 
is, for the bid-ask spread: 

This spread can be seen as a measure of the effect of market incompleteness due to jump risk on the 
price of a specific option, from the point of view of a specific market agent. It decomposes into a 
product of three factors, each representing a specific feature of our market model: 
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Figure 2: Left: jump risk sensitivity as function of K , Sq = 1, a = 0.2. Right: jump risk sensitivity 
as function of T, So = 1, a = 0.4. 


• The parameter a , which characterizes the risk aversion of the economic agent; 

2 

• The factor 7714 — -rj which characterizes the specific Levy model through its variance, skewness 
and kurtosis; 

• The expectation of the integral, which characterizes the specific option, and only depends on 
the variance of the price process. 

The factor 


E 


BS 


d 2 Pj 


BS 


dS 2 


(t, S t ) dt 


(15) 


can therefore be seen as a model-independent measure of the sensitivity of a specific European option 
to jump risk, in the limit of small jumps. It is therefore interesting to study the dependence of this 
measure of jump risk sensitivity on strike and time to maturity. 

Figure [2] plots the expectation (151 as function of strike (on the left graph) and as function 
of time to maturity (on the right graph) for a European put option. We see that the sensitivity 
to jump risk is maximal for options close to the money, since for far from the money options the 
exercise probability and therefore also the price and the spread are very small (remember that we 
are interested in sensitivity to small jumps). Note that actual bid-ask spreads in option markets 
exhibit similar patterns with a maximum close to the money (see Figure [3]), although of course 
actual bid-ask spreads are influenced by a multitude of factors other than jump risk. 

In terms of time to maturity for options which are not at the money, the sensitivity first grows 
(because the exercise probability increases) and then decays for large maturities due to a ‘central 
limit theorem’ effect which smoothes out the effect of jumps. 
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Figure 3: Bid-ask spreads of options on S&P 500 index, observed on Jan 21st, 2006, as function of 
strike, for time to maturity T = 42 days. The underlying index value was So = 1261.49. 

A Proof of Lemma [3] 


Let P\(t,S ) = E* [h(S£(X x )T-t)}- By Proposition 2 in [10], under the assumptions of Theorem 
[2} P\ (t, S) is infinitely differentiable in t and in S, and the assumptions of the Lemma imply, in 
particular, that < L a.s. for all t € [0, T]. Using the martingale representation of the option 

price given in m, we obtain that 


Mji = 


* dS 


crS+dWt 


{zd?S?_ - P x (t , Sl{ 1 + 2)) + P x (t, S*_)}J x .(dt dz ) 


and that the quadratic hedging strategy is given by 

“ = ~dS 




a2 + I R z2iy \(dz) 




zE^{z)v x (dz), 


where we denote 


-—’A 


t (*) = Px(t, 5 t A (l + z)) - P x (t , St) - zS> 


dP x 
dS ‘ 


By the Burkholder-Davis-Gundy inequality expressed in predictable terms m, for q > 2, there 
exist c q ,C q > 0 such that: 


Mi 


E[( sup \M t \Y)<C q E[{M)l + \x\ q *v¥ 

0 <t<T 


(16) 


where is the compensator of the jump measure of the process M. The quantities appearing in 
the right-hand side are explicitly given by 


(M x )t — 


\x\ q *v™ = 


a 2 {Stydt + 


dS 


{zd^St - P x (t , St_( 1 + z)) + P x (t, St_)} 2 v x (dz)dt , 


fj 


Z0is? - P x (t, St_( 1 + z)) + P x (t, St_ )I q v x {dz)dt. 
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Substituting the expression for d x into the first of the above equalities, we further obtain: 

(M a ) t = |j^ dt ^zE x (z)u x (dz)\ + J^zE x {z)v x {dz)-E x {z) S j v x (dz)dt 

rp \ 2 'jp r JT' 

= -^2 f dt \ [ z ^t ( z )v\(dz)j + [ [ (E x (z)) 2 v x {dz)dt < [ [ (Z x (z)) 2 v x (dz)dt. 
a Jo \J r / Jo J r Jo Jr 

(17) 

Our first goal is to estimate ( M x )t■ To this end, we fix 77 e (0, T) and estimate separately 
(M x )T-n and (M x )t — ( M x )t- v ■ 

We proceed with the estimation of (AI x )t- v - By Taylor-Lagrange expansion 

r 1 2 p 

~ x {z) = z 2 (S x ) 2 y o -^(t, P t A (l + 0z))(l - 0)M. 

Note that Vz > —1, j^§^d9 = ( 1+z ) lo g( 1+z l~ z < \ Thus denoting T x (t, S ) = ^Jy-(i, S) 

|S A (z)| < z 2 S x max. S\T x {t, S')|. 


We conclude that 


(M x ) T -n<\ 2 f z A u(dz) f V (S x ) 2 (max S\T x (t, S)\Y dt. (18) 

Jr Jo ^ 

Furthermore, the term ( AI x )t — ( AI x )t- v is controlled by the option’s delta. Indeed, as 


| S | < L, we deduce that 

|P A (t,5 A (l + Z ))-P A (t,S A )| = 

so that finally 

(A1 x ) t - (M x ) T _ ri < [ T [ ( E x (z )) 2 v x {dz)dt 
JT—ri Jr 


J~ S^(t,s x (l + u))du 


< 


m 

1 -S’ 


<2 f f \P X (t, S x (l + z)) — P x {t, S x )\ 2 v x (dz)dt + 2 f [ z 2 (S x ) 2 ° Px 
Jt~ r) Jr Jt—t/Jr 


v x (dz)dt 


(19) 


Using (18 1 , (191 and the gamma estimate of Lemma [ 6 ] we conclude that there exist two constants 
A and B such that 


(TV J x )t P 


X 2 A sup 0 < t < T (5 A ) 2 (logT - log 77 ) + 77 P 
A 2 Tsup 0 < t < T (P i A ) 2 (log T — log 77 ) + - ^ 


TP' 


a > 0 , 
+ 77 P a = 0. 
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Taking r\ = A 2 and assuming A < \/T A b, yields, for a constant C < oo, 


(M x )t < A 2 C \ 1 - sup (5^) 2 log A L 
L 0 <t<T J 

Since S x is a martingale, by the BDG inequality we have, for all q > 2, 


<E[ sup I S x -s x \ q ] 

0 <t<T 


< C q E 


a z + / afv 


(dx)^ ^J S 2 dt^j + A 9 2 J z q v(dz) J S q dt 


which means that E[sup 0<t<T l^l 9 ] is bounded uniformly on A, and therefore, 


E[(M x ) q T /2 ]=o(\ q (logi 


Lastly, it remains to estimate the second term in the BGD inequality (16). In a similar manner 
as for the quadratic variation of M x , we obtain for every q > 2 and for some constant C < oo (which 
may change from line to line), 




xr * v, 


I' I ~ 2 f zE x (z)Mdz)-E x (z) 
Jo Jr a Jr 


v\(dz)dt 


<2 q V 2q \ q 2 f \z\ q v{dz) f f z^ x (z)iz\(dz) dt + 2 q 1 I f \z> x (z)\‘ v\{dz)dt 
Jr Jo Jr Jo Jr 


< CX 2q ~ 2 


r T- v 


| S x maxS\T x (t, S)\\ q dt + Cry X q 


-2 


< 


A 2 <3- 2 c sup 0<t<T (5 t A ) q -hr) +dA q ~ 2 C 

\V 2 T 2 j 


y 2 q~ 2 c SUp 0<t<T (S) A ) q 


% ... . , J 1 + \ 2 (f -1 ) ( 1 . 1 

[-1 T f-1/ +A l 


cr > 0, 


+ qX q ~ 2 C a = 0. 


Choosing once again r) = A 2 leads to, 

\x\ q * i/" A ' p < CA 9 (1 + sup ( S x ) q ) 

0 <t<T 

for some constant (7, and therefore, 

E[\x\ q *u™ X ’ P ] = 0(\ q ) 
as A —> 0. The proof of Lemma [3] is complete. 

Lemma 6 (Estimation of the gamma). 

• Let a > 0. Then there exists C < oo such that 

maxS\T x (t,S)\ < ~^==- 
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> 0 for some /3 G (0,2). Then there exists C < oo such that 


• Let cr = 0 and lim inf 

rj.0 


f[- r , r ] x 2 u{dx) 

7”2 — /3 


maxS\T x {t,S)\ < C ■ 




-l 


y/T-t (T-t)? 


Proof. Let 

Calif (t, S) = E[(S£(X A ) T - t - K) + }. 
Under the assumptions of Theorem [2] 


h{S) = h(0) + ti(0)S+ f ( S-K) + n(dI <), 

J ('0,ool 


(0,c 

where /r is a finite measure on R + , defined by /-/.((a, 6]) = h'(b) — h'(a). Since Calif (f, S) < S, by 
Fubini’s theorem, 

P x (t 1 S) = h(0) + ti(0)S + [ Calif (f, S)ii(dK). 

J (0,oo] 

Moreover, since is bounded and /r is a finite measure, the dominated convergence theorem 

yields 


dP\(t, S) 
dS 


= h'{ 0) 


' (0,oo] 


dCallf (t, S) 
~dS 


n(dK). 


Using the Fourier transform representation for call option price in exponential Levy models [33], 
we get the following identity: 


d 2 Callf ± 
dS 2 1 } 2 tt 


K iu+1 S~ iu - 2 ^_ t (~u)du 

K iu S -iu-i^\_^_ u _ ^ du 


where $ A (u) = E[ e iul °g£(^ A )t] = e *M«) and if x {u) = -^(u 2 +iu) + f (■e™ 1 °gU+*) - 1 - iuz)v x {dz). 
Therefore, 


S ) 

as 2 • ’ b> 


- I ,-<>!*•, 


and the dominated convergence theorem yields 

S\T x (t,S)\<^~ I \^_ t {- u -i)\du=^- f e( r -W*(«-i) du 

Jr 2tt j r 

where C = f (Q , \dfi\ and fif) X {u — i) = + /{(l + x) (cos (u log(1 + a;)) — 1 )}u x {dx). 

Let us study separately the cases a > 0 and a = 0. When cr > 0, we directly get 


$hp x (u — ?)<-- 
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which leads to 


S\T x {t,S)\ < 


C 


When cr = 0, using (201, we get, 


C 


s|r A (t,s)|< — 

\J\u\<t/\ 


< C ( / e-« T -^ 2 du+ I e-imP-Vdu 


cttJ2tt(T — t) 
e -c (T-t)u' du+ f 


o \2-f}( T ([ U 


2tt 

= C 


\u\>i/\ 

ili 


2 1 

1 A? -1 


VT-t (T - t) 


1 ( 5 

P 


where C is a constant. 


□ 


J* is(dx') 

Lemma 7. Assume that a = 0 and liminf 1 ~ r ’ rl 2 _g - > 0 for some /3 £ (0, 2). Then, for every 


rj-O 


C — cu 2 |nA| < £ 
C-c\u\P\P ~ 2 \u\\>£. 


( 20 ) 


£ > 0 and d£R, there exist c > 0 and C < oo such that for u £ 

$bp\(u + iu) < 

Proof. Observe first that 

+ iv) = A” 2 J ((1 + \x)~ v cos(u(l + Aa;)) — 1 + v\x)v(dx) 

= \~ 2 J (l + \x)~ v (cos(u(l + \x)) — l)is(dx) + v(v + 1) J x 2 v{dx) J (1 + d\x)~ v ^ 2 (l — 6)dd. 

Since X has jumps bounded by S and A < 1, there exists C < oo and c > 0 (depending on v) such 
that for all x in the support of v , 

|u(u + 1) J x 2 v{dx) J (1 + 0\x)~ v ~ 2 (l — 0)d6\ < C and |1 + Xx\~ v > c. 


Then, using 1 — cos(:r) = 2(sin |) 2 > 2(^) 2 for |i| < tt, we get: 


5R'0a(w + iv) < C — c [ 
Jh 


n 2 (log(l + Ax)) s 


\U log(l + Atc)|<7T 


A 2 


v(dx), 


but since (log(l + a:)) 2 > £ 2 (log2) 2 for \x\ < 1 we have, for a different c > 0, 


+ iv) < C — c / (ux) 2 l\\ x \<\y{dx). 

J|u log(l+A:r)|<7r 


Once again, by the bound on the jumps of A', | log(l + x)| < on the support of v and we 

also have, for a different c > 0, 


$tip\(u + iv) < C — c j 

J\x\< 


■’(dx). 


AM log(l + «) 


A1 
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J' x^ is (dx') 

Under the assumption lirninf l ~ r ' r] 2 _ fl - > 0 for some /? £ (0,2), there exist rg > 0 and Cq > 0 

r4-0 r 

such that for all r < 7*o, _/|_ T . r .j x 2 v{dx) > cor 2 ~' 3 . This implies that one can find constants £ > 0 
and yet another c > 0 such that 


%hl>\ (u + iv) < 


C- 


i\\ < i 


C-c\u\P\P- 2 \u\\>i 

Now, by changing c and C this inequality can be shown to be true for arbitrary l. 


□ 


B Proof of Lemma H] 

Using the notation of the proof of Lemma [3j we have 


E 


- h(S$) 


= E[(M£) 2 ] = E[(M a )t], 


where ( M x )t was computed in ©■ From the Fourier transform formula for the call option price 

Kr N 1 [ K iu+1 - R Sr iu+R $* J-u-iR) 

Callf = — / - *. , - -du, 

Z7T ,/ R (K — IU)(R — 1 — 


we deduce that 


Calif (t, 5(1 + 2 )) - Call)) (f, 5) - zS 


,5Call f(i,5) 


K/ 


dS 


-S /* 


iii+1 —i? o—iu+ii^A 


>l ) T _ t (—u — iR ) 


. (1 + - 1 + (m - ij)z 


(f? — iu)(R — 1 — m) 


du. 


Since the fraction under the integral sign is analytic for z > — 1, we can choose i? = 1 in this formula, 
obtaining 


Calif (i, 5(1 + z)) - Call)) (f, 5) - z5 


K/ 


,dC all A (f, 5) 
dS 


= — [ K iu Sr iu+1 $T- t (-u - i) ^ + ^ >U+ ! ———*». 

27t J r m(ra - 1) 


and therefore, 


Z?(z) = j pidK) f K iu Sr iu+1 K-t(~u - i) (1 + 2:) ~ ,U+1 . ^ ^ l ^ z du . 

2^ J (0,oo) JR lU(lU 1) 

By Fubini’s theorem, the expression under the time integral in © equals 

At [ n(dK) f n(dK) [ dv [ duK iu K iv 

4tt 2 7(0,oo) 7(0,oo) 7 r 7r 

$T-t(- M -*)4 , T-t(-' y -*)^(-M-v-2i) r a A (u)a A (u) 


—uv(iu — l)(iv — 1) 


+ b x (u,v) 
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where 


a\(u) = / v\(dz)z{(l + z) lu+ — 1 + (iu — l)z} = A / z/(cte) 2 :{(l + \z) lu + — 1 + (iu — l)Xz} 
Jr Jr 

and b x (u, v)= f v\{dz){(l + z)~ iu+1 - 1 + (iu - !)«}{(! + 4 _i,,+1 - 1 + (tv - 1)4 


= A - / j/(d*){(l + A4“ ra+1 - 1 + (iu - 1)A*}{(1 + Xz)~' lv+ - 1 + {iv - 1)A z} 

Jr 

From the explicit form of <h A , 

2 r 

$ t A (u) = e tMu \ ipx(u) = ---{u 2 + iu) + A" 2 / ((1 + A z) iu - 1 - iuXz)v{dz) 

2 Jr 

= — C -—(u 2 + iu)) — (« 2 + iu) f z 2 u(dz) f ((1 + Xz9) vu ~ 2 (l — 9)d9, 
2 Jr Jo 

we deduce, using the fact that v has bounded support, that for every u £ C, 


ip x {u) -> ip°{u) = - 


u z + iu f 2 


+ / z 2 u(dz) I as A-> 0. 


On the other hand, since 


(1 + Xz) lu+1 — 1 + (iu — l)Xz = iu(iu — l)X 2 z 2 f (1 + A 9z) lu 1 (l — 9)d9, 

Jo 


we get that 
A -1 cia(w) 


-i„ x iu(iu- 1) f 3 ,- 2 . v uv(iu- l){iv- 1) /' 4 


z i v(dz) and A ~b\(u, v) —> — - 


z 4 p(dz) 


as A —> 0. Therefore, provided that we can find an integrable bound to apply the dominated 
convergence theorem, X~ 2 (M x )t converges to 


TO4 


m 2 \ 1 


a 2 ) 4-7T 2 


' (0,oo) 


li(dK) / n(dK) / dt dv duK iu K l 


' (0,oo) 


x u — — i)$t(—u — v — 2i) 


= T m 4 it E 


3 \ W BS 


rT 4 2 dP BS {t,S t ) 
5 ‘ dS 2 


dt 


as A y 0. 

We first consider the case a > 0. Remark first that by the bound on the jumps of X, 


[ (1 + A 9z)~ iu£ * 1 {l-9)d9 
Jo 


< 


1 


2(1-4’ 


so that it remains to find an integrable (in u, v and t) bound for 

- i)$th u -v- 2 *) 
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However, in this case, 

|$£_ t (-u - i)&$.-t(~v - - v)\ < e -h(T-t)° 2 (u 2 +v 2 )-ht* 2 (u+v) 2 +t*^ 

which is integrable since 

o -lT* 2 v 2 (l-£) _ 2?r 


du / dve -W-t)a 2 {u 2 W)-hta 2 {u+v) 2 = / 


du < 


aWT 2 -t 2 ’ 


Let us now consider the case a = 0. We shall use the bound (20). In addition, 
\b\(u, u)| < A -2 ( / v{dz)\(l + \z)~ m+1 — 1 + (iu — l)\z\ 


x / v(dz )|(1 + A z)~ lv+1 -l + (iv- l)Xz\ 


and it is easy to show, using arguments similar to those used to prove the bound (201 that for some 
constant C < oo, 

J u(dz )|(1 + Xz) lu+1 — 1 + (iu — l)Az| 2 < CX 4 u 2 (u 2 + l)l|Au|<^ + C'A 2 u 2 l|; lll | > ^, 


where the constant t may be taken the same as in the bound ([20]) . Similarly, by the Cauchy-Schwarz 
inequality, 


|aA(w)«A(^)l < c 2 A 


-iu -\-1 


— 1 + (iu — l)Az|‘ 


i/(dz) |(1 + Xz) 


x v(dz)\(l + Xz) m+1 — 1 + (iv — l)Az| 2 
so that to complete the proof it suffices to study the integral 

[ dv ( du v-2i)\ ( l \ Xu < i \ + —^t==)( 1 |a«<z| + 7 ^ 7 ===)- 

Jr Jr Ayl + u 2 Xyt + v 2 

( 21 ) 

We shall decompose it into four terms corresponding to integrals over non-disjoint sets (whose union is 
R 2 ) {|«A| < 2£, |uA| < 2£j, {|uA| > 2 £, |uA| < £}, {|uA| < £, |t)A| > 21} and {|uA| > £, |uA| > £}, and 
show that on the first set one can apply the dominated convergence theorem, and the contribution 
of the three other sets to the limit is zero. On the first set, the integrand is bounded as follows: 

[*£_*(-« - i)*T-t(~v - -v-2i)\ 1 { |„a|</,|vA|</} < e -< T - t ^ +v ^- ct ^ 2 , 

which is integrable in u, v and t. Hence, by the dominated convergence theorem, 


^4 f KdK) f n(dK) 

4?r •/( 0,oo) J (0,oo) 


dt dv duK lu K v 

J\v\\<2£ J\u\\<2t 


<Z>T-t(~ u - i)®T-t(- v ~ -v-2i) r a\(u)a\(v) 


—> — ( m 4 Z7T ) E 


?bs 


— uv(iu — 1 )(iv — 1) 

: 9P BS (t, S t ) 


+ b\(u,v) 


Si 


dS 2 


dt 
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as A —> 0. It remains to show that the other three sets give a zero contribution to the limit. 

On the set {|uA| > 2£, |uA| < £}, the integrand in (211 is bounded by: 

^ _ —c{T—t)\P- 2 \u\^—c(T—t)v 2 — ct\P-' 1 \u+v\P ✓ ^ _ —c{T—t)\P- 2 \u\P—c{T—t)v 2 —ct\^- 2 \v\P 


vT + u 


a/T + u 2 


A_ -c(T-t)X fi - 2 \u\ fl -cT(i l3 - 2 Al)\v\ 2 


VT-\ 

On the other hand, 


du 


'|uA|>2 £ Aa/1 + U 2 


o - c (T-t) A' 3 - 2 M' 3 < 2 f°° ^g-c(T-t) A- 2 | M |' 3 = 1 j- f T ~t \ 

J2t ^ u A /T — t \ A 2 / 


where 


pOO 

m = 2 v^ / 

J ue 1 


du 


-c|«y 


u 


Note that / is a bounded positive function and f(0) —>■ 0 as 6 —> oo. Therefore, by the dominated 
convergence theorem, 


f dt f du f u — i)$T- t (~ v — *)<!>*(—u — v — 2z)| 

Jo J|«A|>2£ J|dA|<^ 


1 


A a/ 1 + : 


0 


as A —y 0. 

The set {|uA| < £, |i>A| > 2 £} can be dealt with in the same manner. Finally, on the set 
{|uA| > £, |uA| > £}, the integrand in (211 is bounded by 

A" 2 


\/(l + w 2 )(l + V 2 ) 


-c(T-t)X p 2 \u\P — c(T—t)\P 2 \v\ fl -ct\u+v\ 2 -, 

° - L A|'u+t>|<£ 


_|_ g-c(T-t)A' 3 2 \u\ l3 -c(T-t)X 13 2 \v\ l3 -ctX 13 2 \u+v\P ^ (22) 


With a change of variable u = v = ^-4^, using the convexity inequality 


x + y 
2 


+ 


x-y 

2 


> cpQxf + \yf), 


the integral of the first term above satisfies 
r r A” 2 


du 


dv 


|uA|>* J\vX\>Z a/( 1 + u2 )(l + V 2 ) 


-^T-DX^uf-ciT-^-^-ctlu+vl 2 -. 

° X\u-\-v\<£ 


< 


< 


dx 


dy- 


A- 


' A|x+y| >2^ JX\x—y\>2t (1 + |x + t/| ) (1 + \x — y \ ) 

r r A" 2 


D -cc /3 (T-i)A /3 2 \x\ f> -cc ft {T-t)X 13 2 \y\ p -ct\x\ 2 


l A|s|<^ 


dy 


dx 


-c{ Cfi Al)T\x\ 2 < 


X\y\>t JX\x\>£ (1 + \x + J/|)(l + \x - y\) 


<^[ dxe~ c ^^ T l x l 2 , 


for some constant C < oo, where we used Lemma 2 in |5j. It is clear that this expression converges 
to 0 as A —^ 0. 


27 






























Finally, using the same change of variable as above, and once again, Lemma 2 in [5], the integral 
of the second term in (221 satisfies, 


du 


dv 


!\u 


< 


X\>t J\v\\>£ \/(l + U 2 )( 1 + V 2 ) 


-c(T-t)\ ft ~ 2 \u\ fl -c(T-t)\ fi - 2 \v\ tl -ct\ fi - 2 \u+v\ 


1-2 


dx 


dy- 


0 -c(c^A1)TA' 3 - 2 |x|' 3 


l\\x+y\>2£ J \\x—y\>2£ (1 + \x + 2/1) (1 + \x — y\) 

<~ [ dxe-< c ^ TXf, ~ 2 ^ =°, f dxe~< c ^ TX ~ 2 ^ =C\i ~ 3 [ 2 dxe~ c ^ x1 ^ 
^ J \\x\>i ^ J\x\>£ J\x\>£\ 


which clearly goes to zero as A —> 0. 

To perform the computation for the put option pay-off, recall that in the Black-Scholes model, 
S t = Soe~ !Z x t+aWt and d FU gg»' t ’ S '* = < K a: ) = T/W' Therefore, 


E 


BS 




K 


2 fT 


2ira 2 J o 

K 2 ' T 


E 


BS 


' l°g(^t/ K) 

, &VT=t 2 , 


dt 

T-t 


27rcr 2 


E 


BS 


/ i°s(.s 0 /K)-e^+vW t 


7 s/T — t 


dt 

T-t 


It remains to perform the explicit integration with the Gaussian density to get the result. 


C Proof of Lemma [5] 

In this proof, we denote X BS = dW t , where IT is a standard Brownian motion independent from 
A' a . For A > 0, define 

f x (t,x) = E [h(x£(X x ) t £(X BS ) T - t )], 

and for e > 0, let 

h E {x) = E[h(x£(X BS ) e )], P e BS (t , x) = E[h e (x£(X BS )T~t)} = h e+T -\x), 
m,x) = E [h e (x£(X x ) t £(X BS ) T -t)]. 


These functions are well defined because X x has bounded jumps and therefore all moments of 
£(X x ) t are finite. Without loss of generality we shall also take So = 1 below. 

From Ito formula, using item [T] of Lemma [8j 


p bs( p , St) — Tbs(0, 1) + 


dP\ 


BS 


dt 


dt - 


dP\ 


BS qX 


dS 


SyadWt 


1 2 / oX \2 S~PbS J, 

r {S ' 1 ~as^ dt 


dPl 


+ E{ P bs(^ + *)) - pe Bs(t, S x _) - zS x ^(t, S x )}J x (dsdz). 
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Taking the expectation and using the Black-Scholes equation and Fubini’s theorem justified by item 
[3] of Lemma [8j we get: 


E[/i e (^)]-Pg s (0,l) 
= e / t 


v\(dz)dt 


P S BS (t, #(1 + z)) - P S BS (t, S X ) - zS - I.2 (5 A)2^5 
= f [\fx(t,l + z)-f\(t,l)-z{d x fx)(t,l)~ Y{dZfx){t,l)\i'x{dz)dt 

J 0 J M v ) 

= L Jr + j Q yl 1 - 6 ') 3 ( 9 x/a)(^ 1 + ^v\(dz)dt 

Define f x (t,x) = (d 3 f x )(t,x) = E[h e (x£(X x ) t £(X BS )T-t )], where h e {x) = x 3 (d 3 h e )(x). Then, 
using again a Taylor-Lagrange expansion, we get for all t G [0, T\ 


/a(*.1) = /a(0,1)- 


d 3 P\ 


BS 


dS 3 


( 0 , 1 ) 


{/a(s, 1 + z) - /a(+ 1) - z{d x fl){s, 1) - y(^/a)(s, 1 )}v e \(dz)ds 
1 A^ 3 

-^-(1 - 8) 2 (d 3 J e x )(s, 1 + Az9)d0v(dz)ds. 


Substituting this representation into the above formula, we obtain 

A m 3 T d 3 P BS 


E[/+(S£)]=Pg s (0,l) 


6 dS 3 


(0,1) +A 2 [ dt I '—"(dz) [ (1 - 9) 3 (d*f x )(t, 1 + Xz6)d9 
Jo Jr d Jo 


pp I dt I ds f '-v(dz) I (1 - 9) 2 (d 3 f x )(s,l + Xz9)d9. 


Note that \/(s,x) G [0,T] x M + 

( d lf\)( s ,x) = 6(d 3 f x )(s, x) + I8x{d$ft){s,x) + 9x 2 (dlfl)(s,x) + x 3 (d%fl)(s,x). 

Now we use item [5] of Lemma [8] to make e go to zero, obtaining 

Ato 3 T d P BS ^ + X 2 f dt j z u ( dz }f (\-0f^f x )[t,l + \ze)dB 
Jo Jr d Jo 


E[h(S$)\ = P BS (.0, 1) + 


6 dS 3 


A 2 m 3 

6 


dt ds l ^—v(dz) / (1 — 0) 2 (9 3 /a)(s, 1 + Az9)d6. 
Jn J tu J ./n 


To finish the proof, we use the dominated convergence theorem (justified by items [3] and [4] of Lemma 
[8]) to show that 

/ dt. [ '—v(dz) [ (l-9) 3 (djf\){t, l + \z9)d9-+ [ dt f ^—^{dz) [ (1 - 9) 3 (d*f 0 )(t, l)d0 

Jo J r 6 Jo Jo Jr 6 J Q 


m 4 r ,q4i vn U TO 4 T d 4 P B S 


24 


-(^/ o )(0,l) = 


24 dS 4 


( 0 , 1 ) 
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and similarly 


[ dt [ ds f ~—v(dz) f (1 — 9) 2 (dlf\)(s, 1 + Xz6)d6 —»■ —(c^/o)(0,1) 
Jo Jo Jr z Jo iZ 

= ^ { 6P bs(1) + 18^(1) + 94 5 J(1) + P^sO)} . 

Lemma 8. Let the assumptions of Lemma^hold true. Then, for all 0 < k < 6 

1. P E BS eC^([ 0,T]x(0,oo)). 

2. ( d k h e )(x ) exists, is continuous and has polynomial growth in x for all e > 0. 

3. {d k fff)(t,x) exists, is continuous in x and satisfies 

\(d k x fl)(t,x)\<C(l + \x\ n ) 


for some n > 0 and a constant C which does not depend on t, e or X. 
4- For all t,x, ( d k f x )(t,x ) — > ( d k h T )(x ) as A — >■ 0. 

5. For all t, x, A, {d k f{){t,x) ->• {d k f x ){t, x) as e —)• 0. 

Proof. Under our assumptions, the random variable 

log (S(X x ) t S(X BS ) T _ t+E ) 

admits a density p x (t,x) which can be recovered via Fourier inversion: 


p e x (t,x) = J- [ e lux ^ s _ t+E [u)^{u)du. 

Jr 


By the bound (201 and the explicit form of & BS , we conclude that the derivatives of p x (t,x ) with 
respect to x of any order are continuous and given by 


d*Pi(t,x) = — J (-iu)*e-™*?i t+ M*n«)du. 

By Jensen’s inequality and Plancherel’s theorem, for any p > 0, 


[ \d k fpl{t,x)\e^dx= [ \d*pUt,x)e^W\e-Wdx< ( [ \d k x p E x (t,x)\ 2 e^+ 2 ^) 2 

Jr Jr \Jr J 

< ^jjd k pl(t,x)\ 2 e^ +2 ^d x y + (^Jjd k pl(t,x)\ 2 e-^+ 2 ^dxy 

1 

= fj v ~ i( JP + 1 )| 2fe | ®T-t+e(v - i(p + 1 - i{p + l))| 2 d^ 

+ j B \ v + *(p + l)| 2fc |$f-t+ e (« + i(p + 0)®t(v + i(p + !))| 2 ^ < °o. 
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Consider for example, the second term. Using the bound (201, it satisfies, for some constant C < oo, 


/ |i; + i(p + l)| 2fc |$ff t+e {v + i{p + l))$t (v + i(p + l))\ 2 dv 
JR 

< (7 /* (1 _|_ \ v \ 2k ') e -(. T ~t+£)° 2v2 - c t\v\ 2 l\ v x\<e-ct\v\ l3 \ l3 - 2 l\ vX \> ec i v 

Jr 

<C [ {l + \v\ 2k )e~^ T - t) ^ v \e- ct ^ 2 + e~ ct ^)dv, 

Jr 


(23) 


which is easily seen to be bounded uniformly on t. Therefore, f R \d k p x (t, x)\e p ' x 'dx is bounded 
uniformly on f, £ and A. 

This means that the function f x is given by 


[ dz h(xe z )px(t, z). 
Jr 


Instead of the function f x we shall, for notational convenience, study the function f x (t, x ) = /|(f, e x ), 
which is therefore given by 

[ dzh(e z+x )p E x (t,z) = f dzh(e z )p e x (t,z - x) 

Jr Jr 

By dominated convergence, using the above estimate, we then get that 

d k Jl{t,x) = {-l) k [ dzh{e z )d k p x {t,z — x) = (—l) fc [ dzh{e x+z )d k x p\{t,z) 

Jr Jr 

exists, is continuous and has exponential growth in x, which means that d k f x (t,x) has polynomial 
growth uniformly on t, £ and A. This finishes the proof of item[3j 

To study the convergence in A, remark that from the polynomial growth of h , \h{e x )\ < e p ^. 
Then, proceeding similarly to the above, we have 

\d k h(t,x)-d k f Q {t,x)\ < C [ e p|z| | d k p x {t, z) - d k p 0 (t, z)\ 


<C( / (1 + \v\ 2k )\^ s _ t (v ~ i(p + l))(*?(v - i(p + 1)) - > - *(P + l))l 2 dv 


BS( 


+ c / (1 + \v\ 2 k )\^s t{v + i{p + l))($*(u + i(p + 1)) - + i(p + l)))\ 2 dv 

\Jr 

Consider for example the second term. It satisfies 

[ (1 + \v\ 2k )\^ s _ t (v + i(p + imt(v + l(p + 1 )) - + i{p + 1 )))| 2 dv 

Jr 

<C [ (l + |u| 2 fc )e-( T - t ^ 2?,2 |^(u + i(p+l))-$f s (u + *(p+l)))| 2 du 

Jr 

Since + i{p + 1)) —> $f s (v + i(p + 1)) for all v as A —» 0, and an integrable bound can be found 
similarly to ( |23| ), we conclude using the dominated convergence theorem that the above expression 
converges to 0 as A —> 0. This finishes the proof of item [4j Other items can be proved in a similar 
manner. 

□ 
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D A general formula for high-order Black-Scholes greeks 

The risk neutral price process dynamics of the stock in the Black-Scholes model with zero interest 
rate and volatility a > 0 reads 

dS t = aSt.dWt, 

where W denotes a standard Brownian motion. In this model, the pricing function of a European 
call option with strike K > 0 and maturity T > 0 given by P : [0, T] x R + —> R + satisfies 

P(t, s) = E[(5t - K)+\S t = s] = s$(<5i(f,s)) - K$(5 2 {t,s)), 

where <f> and ip are, respectively, the cumulative distribution function and the density of the standard 
normal distribution, and the coefficients di and 82 are defined by 

S n ^ l °Sj( + ^-(T ~t) s u , su , /?F — : 

5i(t,s) = - —- -, o 2 (t,s) = di(t,s) - aVT - t. 


cr\/T — t 


We set for n € N, 


d n P 

d n (t, s) = s n -^(t,s). 


The first two cash greeks can be computed by direct differentiation: 

di(t,s) = s$( 8 i(t,s)) 
2 <p(S 1 (t,s)) 


d 2 (t,s) = s 


saVT — t 


For higher order derivatives of European call/put option prices, the following recurrence relation 
holds for all n > 0 : 


d 3 +n (t, s) = '^2 Ck D n -k{t, s)d 2 +k (t , s) 


k —0 


where C% = (^) = fc! , are the binomial coefficients, 


Dk(t,s) = (—l) fc+i fc! 




and 5(t, s) = 


8 i(t,s) 
o\jT — t 


1 . 


cr 2 (T — t) J 

This recurrence relation leads to the following formulae for the cash greeks up to order 6 : 
d 3 (t,s) = - d 2 (t,s)d(t,s ), 

d 4 (t, s ) = d 2 (t, s ) ^S(t, s) - ^ - d 3 (t, s)S(t, s), 

d 5 (t, s ) = -d 2 (t, s ) (28 (t, s) - + 2 d 3 (t, s ) ^S(t, s) - -J—'j - d 4 (t, s)S(t, s), 

d e (t,s) = d 2 (t, s) (& 6 {t,s) - 2 -^ - 3 d 3 (t,s) ( 2 8(t,s) - 


<J T 


+ 3 d 4 (t, s ) 5{t, s )- 5 - - d 5 (t, s) 6 (t, s). 

(J Z T 


(J*T 
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